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ABSTRACT 

Telescopes like Herschel and the Atacama Large Millimeter /sub millimeter 
Array (ALMA) are creating new opportunities to study sources in the far in- 
frared (FIR), a wavelength region dominated by cold dust emission. Probing 
cold dust in active galaxies allows for study of the star formation history of ac- 
tive galactic nuclei (AGN) hosts. The FIR is also an important spectral region for 
observing AGN which are heavily enshrouded by dust, such as Compton thick 
AGN. By using information from deep X-ray surveys and cosmic X-ray back- 
ground synthesis models, we compute Cloudy photoionization simulations which 
are used to predict the spectral energy distribution (SED) of AGN in the FIR. 
Expected differential number counts of AGN and their host galaxies are calcu- 
lated in the Herschel bands. The expected contribution of AGN and their hosts 
to the cosmic infrared background (CIRB) and the infrared luminosity density 
are also computed. Multiple star formation scenarios are investigated using a 
modified blackbody star formation SED. It is found that FIR observations at 
~500 /zm are an excellent tool in determining the star formation history of AGN 
hosts. Additionally, the AGN contribution to the CIRB can be used to determine 
whether star formation in AGN hosts evolves differently than in normal galaxies. 
The contribution of Compton thick AGN to the bright end differential number 
counts and to the bright source infrared luminosity density is a good test of AGN 
evolution models where quasars are triggered by major mergers. 

Subject headings: galaxies: active — galaxies: evolution — galaxies: quasars: 
general — infrared: galaxies — X-rays: galaxies 
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Introduction 



T hanks to the Spitze r Space Telescope and its imaging and spectroscopic instrument 
suite ( IWerner et al.l |2004| ) . it has recently been found that active galactic nuclei (AGN) 
can dominate the near and mid-infrared luminosity an d make a significan t contribution to 



the 8-1000 Atm luminosity, Ltr, of AG N host galaxies (iHouck et al.l 120051 : lYan et al.l 12005 



Weedman et al.ll2006t ISoifer et al.ll2008l . and references therein). Specifically, a large fraction 
of ultra-luminous infrared galaxies {Ljn > 10^^ L©; ULIRGs) and luminous infrared galaxies 
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to determine the AGN contribution to the combined AGN and host observed flux as little 
is known about the AGN spectral energy distribution (SED) at these wavelengths. For the 
purposes of this study we define the FIR as 70-1000 /im. 

The next several years will see a huge growth in high quality FIR and submillime- 
ter spectral data due to telescopes like the Herschel Space Observatory and the Atacama 
Large Millimeter /submillimeter Array (ALMA). Herschel, launched in May 2009, offers un - 
precedented coverage of the FIR and submillimeter spectral regions (jPilbratt et al.ll2010l ). 
Herschel carries two imaging photometers, the Photodetector Array Camera and Spectrorn- 
eter (PACS) which has photometric bands at 70, 100, and 160 /xm (jPoglitsch et al.ll2010l ) 
and the Spectral and Phot ometric Imaging R eceiver (SPIRE) which has photometric bands 
at 250, 350, and 500 yum (iGriffin et al.ll2010f). As one of the p rimary science goals of Her- 
schel is to study the evolution of galaxies (jPilbratt et al.ll2010l ). three wide-area surveys of 
various depths are being undertaken. The PACS evolutionary probe (PEP) will cover ~3 
deg^ and will observe some fields, such as th e GOODS fields, d own to the depth of the 100 
and 160 yum confusion limits of a few mJy (iBerta et al.ll2010l ). The Herschel Multi-tiered 
Extragalactic Survey (HerMES) will o bserve ~70 deg^ with SPIRE, reaching 5c r depths of a 
couple mJy at 250 /im in some fields (lOliver et al.ll2010l : iRoseboom et al.ll2010l ). Observing 
~510 deg2 with both PACS and SPIRE, the Hersc hel Atlas (H- ATLA S) survey is the largest 
open time survey being undertaken by Herschel ( Eales et al. 201oh. H- ATLAS will reach 
5a depths of 67 mJ y at 100 /im and 53 mJy at 500 /im (lEales et al.ll2010l : llbar et al.ll2010l : 
Pascale et al.ll2010l ). The ground-based instrument Submillimeter Common-User Bolometer 
Array 2 (SCUBA-2) will operate at 450 and 850 /im and will off er an intermediate se nsitiv- 
ity and mapping speed between those of Herschel and ALMA (iHoUand et al.ll2006l ). Also 
ground-based, ALMA is an imaging and spectroscopic instrument operating in the millimeter 
and submillimeter regime. The ALMA primei0 notes that at full science operations, starting 



^available at http://www.alinaobservatory.org/iinages/stories/publications 
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in late 2012, ALMA will have a 450 /im band (Band 9) with a continuum sensitivity of 0.69 
mJy and a 350 fim band (Band 10) with a sensitivity of 1.1 mJy, along with bands at longer 
wavelength^. The 450 fim band will be available for early science operations in mid 2011, 
albeit at reduced sensitivity. Herschel^s observing capabilities, and in particular the SPIRE 
deep and wide surveys at 250, 350 and 500 /xm, will provide promising candidate targets for 
ALMA early science operations. 

It is expected tha t AGN and their hosts will be readily detected in the FIR by Herschel. 
Hatziminaoglou et al.l ( l2010l ) found that a third of their AGN sample had Herschel 5a de- 
tections at 250 /im with f25o^im > 12.8 mJy, where f250fj.m is the 250 fim flux. The FIR flux 
from AGN and their hosts will be due to dust heated by a combination of AGN radiation 
and star formation. Observations suggest that only 3-9% of subm illimeter galaxies are dom- 
inated by AGN emission in the submillimeter (ILaird et al.ll2010l ). Howe ver, a large portio n 
of galaxies, nearly 50% of the galaxies in the infrared selected sample of lYuan et al.l (120101 ). 
show strong infrared emission from both star formation and AGN activity. Consequently, 
the AGN contribution to the FIR emission of AGN hosts must be considered in order to not 
over-estimate star formation rates in AGN hosts. 

As the FIR SEP of AGN is expected to g enerally be dominated by star formation (e.g. 
Hatziminaoglou et al.l l2010l : iLutz et al.l 120101 ) . FIR observations will inform investigations 
of the star formation history of AGN hosts. By comparing the star formation history of 
AGN hosts and normal galaxies, the role of AGN in galaxy evolution will be constrained. 
Understanding AGN host star formation histor y is therefore an impo rtant tool in determin- 
ing what processes trigger AGN. For example, ISanders et al.l (119881 ) proposed that galaxy 
mergers trigger both intense starbursts and quasar activity. Recent observations and sim- 



ulations support this mode l (IPage et a. 



Draper fc Ballantynd l2010l : iKelly et al 



20041 : iHopkins et a. 



200a iRigopoulou et all 12009 



2010l). IPage et all J2004r found that type 2 AGN 



tend to have higher SFRs than type 1 AGN. jPraper fc Ballantynd (120101 ) found that a sig- 
nificant fraction of Compton thick (CT) AGN, AGN with an X-ray obscuring column density 
N h > 10^^ cm~^, are q uasars which are accreting very rapidly, as predicted by simulations 
by lHopkins et al.l (120061 ) and others. However, it appears that this model may only be appli- 
cable for the most powerful AGN and that Seyfert strength AGN may evolve mor e secularly 
(iBallantyne et al.l l2006al : iHasingerl l2008l : iLutz et al.l l2010l : iNarayanan et al.l |2010| ) . In order 
to determine the applicability of the merger driven evolution and secular evolution models, 
the star formation history of AGN hosts must be accurately determined. This requires a 
robust method of identifying AGN hosts and the AGN contribution to the FIR emission. 



^Sensitivities are for an integration time of 60 seconds and a spectral resolution of 1.0 km/s. 
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FIR emission from AGN comes from the "torus" of dusty gas, which, acc ording to the 



unified scheme, surrounds the central engine of all AGN (e.g., lAntonuccil Il993l ) . This dusty 
torus will absorb higher energy emission from the central engine of the AGN and re-radiate 
this energy in the infrared. The temperature of the dust, and therefore the peak wavelength 
of the infrared radiation, is dependent on the densit y of the obscuring gas an d the distance 



of the obscuring gas from the central engine (e.g., iBallantyne et al. 



2006br). arnong o ther 



parameters hke the geometry of the obscuring torus (e.g., iHatziminaoglou et al.ll2009l ). It 



Elvis et al. 


1994; 


Netzer et al. 


2007) 



reach CT levels of obscuration, it is expected that the clouds of dusty gas in CT AGN will 
have a greater spatial extent than in less obscured AGN. Therefore, CT AGN will have a 
significant reservoir of dusty gas which is cooler than most AGN-related dust clouds. This 
cooler dust component will cause CT AGN to be brighter in the FIR than unobscured AGN. 



Indeed, studies suggest that bright galaxies detected in the F 
heavily obscured, Nh > 10^^ cm~^, and pos si bly CT AGN (e . g.. 



R host a large number of 



Alexander et al. 2005a b 



Bonriovanni et al.lboiol Iwilman et aPboioh . iMuUanev et~al1 (I2OI0I) found that 45-75% of 



the 1 Ms CDF-S X-ray sources will be detected by Herschel at 100 /xm and that deep infrared 
observations with Herschel will allow for a significant fraction of the CT AGN population 
to be identified . Also, since X-ray selection misses nearly half of all infrared identified AGN 
JFu et al.ll20l"ol ) it is expected that FIR telescopes like Herschel and ALMA will detect highly 
obscured AGN that are missed by the deep X-ray surveys. However, ther e is an open debate 
as to how many of these infrared AGN reach CT levels of obscuration. iFiore et al.l (120091 ) 
find that up to 90% of sources with /24 > 550 /iJy and f2i l Jb > 1000 are CT AGN, wher e /24 
is the 24 /im flux and /r is the R-band flux. Conversely, iGeorgantopolous et al.l (I2OIOI ) find 
a more moderat e fraction of CT AGN for this population. When looking at similar infrared 



excess sources. 



Georgakakis et al.l (120101 ) found no strong evidence suggesting that these 
sources host AGN with CT levels of obscuration. Therefore further study of the infrared 
properties of AGN are important in determining exactly how many CT AGN are hidden 
within dusty luminous infrared galaxies. 

This work makes predictions for FIR AGN number counts, contribution to the CIRB, 
and luminosity density using a population synthesis model informed by constraints from the 
cosmic X-ray background (CXRB) and deep X-ray surveys, includin g the Eddington ratio 
dependent CT fraction, fcr, of the composite model investigated by [Draper fc Ballantyne 
(120101 ). The effect of various AGN host star formation scenarios are also investigated. The 
calculation of the model SEDs is discussed in §21 Predictions for bare AGN are presented in 
[|3] followed by predictions for various host star formation scenarios in ||H In §5] and §6] the 
results are discussed and summarized. A ACDM cosmology is assumed as necessary, with 
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ho = 0.7, n„ = 0.3, and = 0.7 flHinshaw et allboOOf ) 



Calculation of AGN SEDs 



As the goal of this study is to make predictions of the average FIR properties of AGN, 
SEDs are used which are representative of an ensemble of AGNs at a given 2-10 keV lumi- 
nosity, Lx, and z, instead of using a SED template based on observations of a statistically 
small set of AGNs. Since AGN IR emission is primarily due to the obscuring gas and dust 
re-radiating absorbed X-ray emission, photoionization simulations allow for the computation 
of the IR emission of an average AGN with a given Lx and Nh- Also, this SED computation 
method provides an opportunity to explore the parameter space of obscuring gas location, 
density distribution, and dust content, which, upon comparison with observations, may offer 
constraints for these physical parameters. 



Similarly to iBallantyne et al.l (j2006bl ) , the IR AGN SED s are calculated using the pho- 
toionization code Cloudy version COS. 00 (IFerland et al.lll998l ). Cloudy includes the compli- 
cated physics of radiative transfer through dusty gas and uses a physical dust model which 
takes into account silicate and graphite grains along with polycyclic hydrocarbons (PAHs). 
This technique does necessit ate a simplification of the IR emitting region, w hich might actu- 



ally be quite complex (e.g., iHatziminaoglou et al.ll2009l : iHonig et al.ll2010l ). However, since 



our purpose is to describe average properties of AGN, and not to model individual objects, 
this technique is an appropriate method of SED calculation. This assumption will be tested 
by comparing the model SEDs against real data in §2.3[ 



2.1. Cloudy Model Setup 



Each Cloudy model is setup such that a constant AGN spectrum, characterized by Lx, 
is incident upon a cloud with inner radius r from the continuum source. The inner radius r is 
set to 10 pc. Compton thin clouds are assigned a uniform hydrogen density uh = 10'^ cm~'^, 
as is typical of molecular clouds. In order to prevent the dust mass from becoming too large, 
the Compton thick clouds are assigned uh = 10^ cm~^. Molecular clouds of this density are 
not uncommon and have been observed in th e Large Magellanic Cloud (IRubio et al.ll2009l ) 
and in the Orion Nebula (jPersson et al.ll2007l ). Simulations were also run with r = 1 pc and 
Uh = 10^ cm~^ for both Compton thin and CT clouds to test the sensitivity of the results 
to these assumptions, which is discussed later. 



The Cloudy model input files are similar to those used by 



Ballantyne et al, 



(j2006b|), with 
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the fo llowing: improvements. Instead of using a constant aox for all AGN, the ISteffen et al. 
( 120061 ) Q:ox-^2keV relation is used here to determine aox- The most significant improvement 
made includes dividing the AGN population into Eddington ratio bins, high (L/Leha > 0-9, 



where L is the bolometric luminosity found using the bolometric correction by lMarconi et al. 



(I2OO4J ) and LEdd is the Eddington luminosity), moderat e (0.9 < L/LEdd < 0.01), and low 
(0.1 < L/LEdd), as described bv iDraper fc Ballantvnd ( I2OIOI ). Both the composite and 
original models of Draper fc Ballantyne ( 201ol ) are investigated. The fax of the composite 
model is Eddington ratio dependent and finds that ~86% of AGN accreting at greater 
than 90% of their Eddington rate and ~60% of AGN accreting at less than 1% of their 
Eddington rate are Compton thick. The original model assumes Compton thick AGN are 
a simple extension of the Compton thin type 2 population and that ~44% of type 2 AGN 
are Compton thick. Since the covering factor is set assuming the unified scheme holds, the 
covering factors are Eddington ratio dependent since for is Eddington ratio dependent. For 
Compton thick objects the covering factor is set as the Compton thick fraction, fcr- For 
objects with 22 < log Nh < 24 the covering factor is se t as (1.0 — fcT)h where fo is the 



type 2 fraction and is calculated as discussed in §2.2 of [Draper fc Ballantynd (120091 ). The 
covering factor for objects with 20 < logN^ < 22 is set to (1.0 — /ct)(1-0 — f^)- As the 
covering factor is dependent on the Eddington ratio and varies from 2; = to 1, Cloudy 
models had to be calculated as a function of z and L/LEdd for each Lx and Nh- 

There is evidence that different levels of obscuration in quasars might be related by an 



evolutionary scenario instead of by orientation effec ts (e.g., iSanders et al.lll988l : iPage et al. 
20041 : braper fc Ballantvnelboiol: bonlev et allboioh . and thus the covering factor need not 
be related to fcT or /s. iBallantvne et all J2006bh found that the unified model assumption 
holds for lower luminosity quasars and Seyfert galaxies but does not seem to hold for high- 
luminosity quasars, possibly due to a different evolution for high luminosity, and therefore 
high Eddington ratio, quasars. Here this issue is addressed by using the Eddington ratio 
dependent fcT of the composite model, which takes into account the different evolution of 
high Eddington ratio quasars from the moderate and low Eddington ratio AGN. 



2.2. The Model Grids 

As this investigation of the FIR properties of AGN is informed by the constraints offered 
by the CXRB and X-ray observations of AGN, SEDs are calculated as a function of logLx 
and z. It is assumed that the redshift evolution of /2 halts at z = 1, thus models are 
computed up to 2; = 1 and the z = 1 models are used at z > 1. Therefore, we compute 
Cloudy models for 2; = to 1, in steps of 0.05, and logLx = 41.5 to 48, in steps of 0.25. 
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For each luminosity and redshift, Cloudy models are calculated for each log Nh = 20.0, 
20.5, . . ., 24.5, 25.0 cm~^. Since the covering factor is Eddington ratio dependent. Cloudy 
models had to be calculated for each luminosity, redshift, and column density for each of 
the three Eddington ratio bins. This resulted in 18711 individual Cloudy models for each r 
investigated and an additional 5103 Cloudy models to investigate the sensitivity of results 
to the hydrogen density, uh, of the Compton thick clouds. 

As we need to run 33 models for each {Lx, z) pair, it is computationally prohibitive to 
run models on a finer grid. Therefore, we have linearly interpolated between the SEDs to 
allow for a finer grid in the calculations. Thus eight SEDs are interpolated between consec- 
utive log Lx model SEDs at a given z and two SEDs are interpolated between consecutive z 
model SEDs at a given Lx- A convergence test was conducted using twice as many steps in 
both Lx and z; it was found that the step size used here is adequate for convergence. 

To create final S EDs fr om the 18711 Cloudy models, we follow the method described 



by iBallantyne et al.l ( l2006bl ). For each Eddington ratio bin, the weighted average of the 
reflection components of each is added to the net transmitted continua of each to 
create the "unified SEDs". The final SEDs are calculated in three categories: type 1, type 
2 but Compton thin (which will be referred to as "type 2"), and Compton thick, for each 
Eddington ratio bin. The type 1 SED is an average of the 20 < log Nh < 22 unified SEDs, 
the type 2 SED is an average of the 22 < log Nh < 24 unified SEDs, and the CT SED is an 
average of the 24 < log Nh unified SEDs. At this point, for each {Lx,z) pair there is a type 
1, type 2, and Compton thick SED for each Eddington ratio bin. Figure [T] shows the rest 
frame logLx = 43 and z = 0.45 SEDs for the high Eddington ratio bin, which has fcT = 
0.86 and /s = 0.78. 



2.3. Properties of AGN Model SEDs 



The model SEDs were tested in the same manner as by IBallantyne et al.l (l2006bl ) to as- 
sure the assumptions used to compute the SEDs were appropriate and that the model SEDs 
are consistent with observed SED trends. As little is known about AGN FIR SEDs, these 
tests primarily investigate the mid-infrared properties of the model SE Ds. The mid- i nfrare d 



colors were examined and found to be consistent with the findings of iBrand et al.l (120061 ). 



(Lutz et al. 


2004; 


Riffbv et al. 


2004 



Gandhi et al.ll2009l ). Also, the fluxes were tested against and found to be consistent with 



expec tations from local bright AGN, known as "Piccinotti AGNs" (see lAlonso-Herrero et al. 
20041 ). As our goal is to make predictions which characterize an ensemble of AGNs, like 
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number counts, the model SEDs are not compared against individual objects but against 
data which explores average properties of AGN SEDs. The model SEDs are found to be 
consistent with the observed AGN SED trends. 

The model SEDs were also tested against the more recently discovered c orrelation be 
tween mid- infrared luminosity, Li2.3^m, and Lx for a sample of local Seyferts (iGandhi et al. 
20091). Figure | 2] com pares the z = 0.05 model SEDs to the local best-fit correlation found by 



Gandhi et al.l ( 120091 ). The model SEDs agree reasonably well with the correlation, especially 
the type 1 SEDs. At high luminosities the Compton thick SEDs appear to have a slight 
X-ray excess compared to the local Seyfert correlation. This is due to the Compton thick 
SEDs including more power radiated in the FIR, and consequently less power radiated in 
the mid-infrared, than the Compton thin SEDs. The r = 1.0 pc model SEDs do match 
the Li2,3fj,m ^ Lx relation better than the r = 10 pc model SEDs; however, the r = 10 pc 
model SEDs best descri be the general mid-infrared properties of an ensemble of AGN (see 
Ballantyne et al.ll2006bl ). The fact that the r = 1.0 pc model SEDs fit this relation better 
than the r = 10 pc model SEDs shows that the distribution of gas and dust in the AGN 
dusty torus is an importa nt effect in the mid-i nfrared. In order for the models used here 
to better approximate the iGandhi et al.l (120091 ) relation, the geometry of the dusty torus 
should be carefully taken into account. The FIR SEDs are dominated by star formation, so 
the geometric details of the hot dust close to the central engine should not affect the results 
of this study. Also, as this study is focused on the integral properties of a large number of 
AGN, detailed modeling of the torus geometry is beyond the scope of this current study and 
we leave this for future work. Overall, the model SEDs are an appropriate representation of 
average AGN SEDs. 



Predictions for Bare AGNs 



In this section we present predictions for differential number counts, the AGN contribu- 
tion to the cosmic infrared background (CIRB), and lun ainosity density, based on th e model 
SEDs. The population synthesis model last described by [Draper fc Ballantynd ( l2010l ) is used 
to incorporate the information about and the constraints on the AGN population from deep 
hard X-ray surveys and the CXRB into infrared predictions. 
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3.1. Differential Number Counts 

The number of sources per square degree witli flux greater tlian 5*, A^(> 5*) is found by 

^^0 Jz^,„ imax(logL--,logLS) d log Lx (1 + z)^ [iZ^ (1 + 2;)^ + JIaJ 

(1) 

wfiere tlie factor Kf^^ = 3.05 x 10~^ converts from s r~^ to deg~^, d^\(Lx, z)/d l og L y is the 
evolving Eddington ratio space density computed by lDraper fc Ballantynd (120101 ). in Mpc~^, 



di is the luminosity distance, and log Lx is the 2-10 keV rest-frame luminosity corresponding 
to the observed-frame infrared flux S at redshift z. The differential number counts are found 
by taking the derivative of A^(> S) with respect to S, dN{> S)/dS. The differential number 
counts were also calculated directly. However, due to the coarseness of the Lx and z grid 
used, the differential number counts contained numerical artifacts at brighter fluxes. The 
predictions from the two calculation methods are in agreement; however, computing the 
integral number counts on a very fine flux grid and then taking the derivative with respect 
to flux minimizes the numerical noise found in direct calculation of the differential number 
counts. 

The Euclidean normalized differential number counts for 70, 100, 160, and 250 /im are 
shown in figure [^. Plots are not shown for the 350 and 500 /xm bands because, as evident in 
figure [H the flux due to the AGN is very small at such long wavelengths; for example, at 10 
mJy the 500 /im band dN{> S) /dS ^ 0.4 AGN hosts mJy^-^ deg~^. The black hues show the 
predictions for the composite model and the cyan lines show the predictions for the original 
model. For ease in interpreting the figure, the low, moderate and high Eddington ratio bins 
have been combined and the total differential counts are shown (see §5.41 for a discussion 
about the contribution from the various Eddington ratio bins). As expected, the differential 
number counts for bare AGN are very small at long wavelengths because the AGNs tend to 
create hot dust. Therefore, long wavelengths can be used to investigate the star formation in 
AGN hosts. For both the original and composite models, CT AGN dominate the bright end 
of the differential counts, especially at wavelengths longer than 70 /im. In both models the 
type 2 AGN dominate at lower fluxes and the type 1 AGN contribute significantly less than 
their obscured counterparts at all wavelengths. The original model and composite model 
differential counts peak at about the same flux for each wavelength and have similar bright 



■^Fluxes are calculated at the filter nominal wavelengths. Using the full filter transmission function to 
calculate the dN/dS provides a result which is within a factor of 1.05 of that obtained when calculating 
the flux at the filter nominal wavelength. Thus the error due to calculating the fluxes at the filter nominal 
wavelength instead of using the full filter profile is negligible compared to the uncertainties in the model and 
the measurements. 
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end slopes. However, the original model peaks above the composite model and declines faster 
on the lower flux level end. 

It is important to consider the dependence of these predictions on the parameters of 
the distance of the inner edge of the cloud from the ionizing source, r, and the uh of the 
CT clouds, nH,cT- As seen in figure IH when r = 1 pc the differential counts are greatly 
reduced at all FIR wavelengths because the dust is hot and radiating more energy in the 
mid-infrared. Contrastingly, when r > 10 pc the differential counts are enhanced due to the 
dust being farther from the illuminating source and therefore cooler. Similarly, when uh^ct 
= 10^ cm~^ the differential counts are slightly enhanced in the FIR due to the greater spatial 
extent of the cloud. However, since the combined AGN and host SED will be dominated by 
star formation at longer wavelengths, and there is likely to be variation among individual 
objects, the exact values for r and uh^ct used will have little effect on the final predictions. 



3.2. CIRB 



The AGN contribution to the CIRB is calculated similarly to population synthesis mod- 
els looking at the CXRB, but exchanging the X-ray spectrum with the infrared spectrum. 
Thus we have 

c_ Z-'"*^^™^ d^,{Lx,z) ^ S,{Lx,z)d^ 

Ho 



^dlogLxdz, (2) 



nogLT- dlogLx {1 + zy[Q^{l + zf + QaY/' 

where S^{Lx,z) is the observed-frame AGN model SED, in Jy sr^^, computed with Lx at 
redshift z. As there is little knowledge about the SED s of AGN in t he FI R, there are very 
few constraints on the AGN contribution to the CIRB. iJauzac et al.l (120101 ) claim that AGN 
contribute < 10% of the CIRB at 2; < 1.5, based on Spitzer observations in the GOODS 
and COSMOS f ields. By extrapo lating from CXRB models to the CIRB using rough AGN 
SED templates, Isilva et al.l koo4 ] predict that bare AGN co ntribute ^0.3% of the C I RB at 
160 /im0. Using a method similar to the one described here, Ballantyne fc Papovich ( 2007 ) 
predict that bare AGN contribute ~0.9% of the CIRB at 160 /xm. Here it is found that both 
the composite and original models predict that bare AGN contribute ~0.9% of the CIRB at 
160 fim. As seen in figure 0, type 2 AGN dominate the CIRB at lower wavelengths and at 
higher wavelengths CT AGN dominate the AGN contribution to the CIRB. At wavelengths 
greater than ~100 fim the original model and the composite model predict very similar 
contributions of AGN to the CIRB. Below ~100 /im the composite model predicts a higher 
contribution to the CIRB by AGNs than the original model. 



*Here we assume the intensity of the CIRB at 160 /im is 12.84 nW m ^ sr ^ (jAltieri et aLll201Q ) 
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3.3. Luminosity Density 



The cosmic infrared luminosity density and its evolution are good indicators of the 
star formation rate density and the evolution of star formation in the universe. However, 
infrared luminosity density is significantly contaminated by AGN, thus making it important 
to understand the contribution of AGN to the infrared luminosity density and how the AGN 
contribution to the infrared luminosity density evolves. The 8-1000 /xm AGN luminosity 
density, piR, is calculated as 



PIRK.Z) 



-'IR- 



I R,min 



dlogL 



dlogLiR, 



IR 



(3) 



where d^\/d log Ljr is the evolving Eddington ratio space density calculated in lDraper fc Ballantyne 
( 120101 ) in terms of the total infrared luminosity, i.e. a total infrared luminosity function for 
a population of AGN with a specific Eddington ratio, which is found using the relation 

d^x{.LiR,z) ^ d<^x{Lx,z) djlogLx) 
d(\ogLiR) d(\ogLx) d(\ogLiR)' 

As seen in figure El both the composite and original models are in agreement with obser- 
vations by iGoto et al.l (l2010al ) at 2; ~ 0.2-0.8. For z > 1 both models do not evolve fast 
enough to agree with the piR and ULIRG AGN infrared luminosity density, P^r^^^ ■ How- 
ever, both models are in fairly good agreement with the high redshift LIRG AGN luminosity 
density, pYr^^ ■ Locally both the original and composite model under-predict t he ptr and 
over-pr edict the pYr^^'^, but are in agreement with the P^r^'^, as measured by iGoto et al. 



(l2010bl ). Due to the r nethod used to separate the AGN and star forming galaxies in the 
samples of Goto et al. ( 2010a| jbl). there is an uncertainty of up to 50% due to the complica- 
tion of gal axies whose infrare d SED is not clearly dominated by AGN activity nor by star 
formation ( Goto et al. 2010bl ). Also, by applying the same method of separating the con- 
tribution of AGN and star forming galaxies to the infrared luminosity density measured by 
Le Floc'h et al. (j2005 ) using Spitzer 24 pm sources in the CDF-S, the local p^r^*^ is ~0.4 
dex higher and in good agreement with local P^r^^^ predicted here. 



4. Accounting for Star Formation 



The infrared SEDs of most AGN are dominated by star formation (e.g.. lFranceschini et al 



20101). Unless the AGN host galaxy can be spatially resolved, it is very difficult to separate 
the AGN emission from the star formation emission. Therefore it is important to consider 
how various AGN star formation scenarios will affect the predictions discussed above. In par- 
ticular, a constant star formation rate is considered, keeping with the unified scheme. A star 
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forma tion scenario consistent with the AGN evolution scenario discussed by ISanders et al. 
( 119881 ) is also analyzed. Finally, star formation scenarios which include evolution of the SFR 
with redshift and AGN Lx are considered. 

To calculate the star formation SED for a given SFR, a modified black bod y spectrum 
charac terized by the dust emissivity, /3, and the dust temperate, T^, is used. The iKennicutt 
fll998h relation, 

SFR = 4.5 X 10-^%jj, (5) 

is us ed to determine Ljr, in erg s~\ for a given SFR, in yr~^, and the Ljji-Td relation 
from kmblard et all (120101 1 . 

Td = To + a\og{LiR/LQ), (6) 



is used to determine T^. Usi ng Herschel observatio ns, lAmblard et al.l ( 120101 ) find Tq = 20.5 K 
and a = 4.4, using /3 = 1.5. lAmblard et al.l (120101 ) define Ljji as the 8-1100 fim luminosity, 
instead of the 8-1000 /im luminosity as done here; however, for our purposes, the difference is 
negligible. The resultin g star format i on SE Ds are in good agreement with the star formation 
templates presented by iRieke et al.l (120091 ) at rest frame wavelengths greater than 50 /zm. 
The star formation SEDs are then added to the model SEDs discussed above to create the 
AGN+SF SED. The calculations discussed in section §3] are repeated using the AGN+SF 
SEDs in place of the AGN model SEDs. Here, the predictions for the composite model are 
presented. A comparison of the predictions of the original and composite models is discussed 
in 



4.1. Constant Star Formation 



The simplest star formation scenario is where AGN hosts have a constant SFR. lBallantyne fc Papovich 



(120071 ) found that an average AGN SFR = 1.0 Mq yr~^ reproduces the AGN contribution 
to the mid-infrared portion of the CIRB as measured by Spitzer, thus for the constant star 
formation model we set SFR = 1.0 Mq yr~^. In keeping with the unified scheme, all AGN 
have the same SFR regardless of spectroscopic type. The Euclidean normalized differential 
number counts for the constant star formation model are shown as the black lines in figures 
El El m Uni [m and At wavelengths shorter than 350 fim, the number counts increase 
steeply from 0.01 mJy to ~1 mJy and then remain approximately flat through ~25 mJy. At 
wavelengths longer than 350 /im, the number counts turn over around 1 mJy and continue 
to decrease toward brighter flux levels. At all wavelengths the CT AGN dominate the bright 
end counts above ~1 mJy, with type 2 AGN dominating at lower flux levels. 



The black lines in figure [13] show the AGN contribution to the CIRB for the constant 
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star formation model. In this star formation scenario, AGN contribute ~4% of the CIRB 
at 160 yum. The AGN contribution to the CIRB is dominated by CT AGN above 100 fim 
and type 2 AGN also contribute significantly. The AGN contribution to the CIRB peaks at 
approximately the same wavelength as the CIRB itself peaks, ~160 fim. 

The pm is shown in figure [TH with the constant star formation model shown as the black 
lines. The p^ji^^^ is not increased over the bare AGN scenario, however, the total and LIRG 
AGN infrared luminosity densities are increased slightly. Thus the constant star formation 
model does not allow for rapid enough evolution with redshift to match the observed high 
redshift piR. 



4.2. AGN Evolution Scenario 



It has been suggested that, at least in the quasar regime, the level of o bscuration ob 



served in an AGN is connected to the evolutionary stage of the qu asar (e.g., ISanders et al. 



19881 : IPage et al.ll2004l : lBallantynell2008l : [Draper fc Ballantynell2010l ). Galaxy merger simula- 
tions support the evolutionary scenario showing that gas rich m ergers lead to a burs t of star 



formation and intense, highly obscured black hole growth (e.g., iHopkins et al.ll2006[ ). In the 
AGN evolution scenario, high Eddington ratio CT AGN would be expected to have SFRs 
reaching into the ULIRG regime and type 2 AGN hosts wou ld be expected to have more 
star formation that type 1 AGN hosts (e.g.. iPage et al.ll2004l ). Following this prescription 
we set the high Eddington ratio SFR = 175 M© yr'^ the type 2 SFR = 2.0 M© yr'^ and 
type 1 SFR = 0.5 Mq yr~^. The low Eddington ratio CT AGN are also given SFR = 0.5 
Mq yr~^ as these are weak AGN mos t likely obscured by mole cular clouds in the host bulge 
or by dust lanes in the host galaxy (IMartmez-Sansigrd 120091 ). and not by intense nuclear 
starbursts. This star formation scenario gives predictions which are only negligibly different 
from the scenario where all AGN have SFR = 1.0 yr~^ except for the high Eddington 
ratio CT AGN which have SFR = 175 Mq yr"!. 

The Euclidean normalized differential counts for the AGN evolution scenario are shown 
as the green lines in figures [TJ [H [9l [TOl [TTl and [121 At all wavelengths the number counts 
rise as the flux level increases, leveling off around 1 mJy, and then continue to rise at least 
until the flux level of 10 mJy. CT AGN dominate at flux levels greater than ~6 mJy and 
type 2 AGN dominate at lower flux levels for all wavelengths. 

In figure [131 the AGN contribution to the CIRB for the AGN evolution scenario is 
shown in green. At 160 yum, AGN contribute ~6% of the CIRB. CT AGN dominate the 
AGN contribution to the CIRB at wavelengths greater than ~200 /xm and type 2 AGN 



- 14 - 



dominate below 200 /iin. The peak of the AGN contribution to the CIRB is roughly at the 
same wavelength as the peak of the CIRB. 

The green lines in figure [T3] show the pm for the AGN evolution scenario. The total 
PiR is in decent agreement with observations except at the highest redshift bin, which this 
model under-predicts. The local P/^^*^ and f^ji^'^ are over-predicted by this model but are 
in good agreement with observations at higher redshifts. However, the P/'r^^*^ observations 
show a stronger redshift evolution than predicted by this model. 



4.3. Evolution with Redshift 



Also considered was t he star formati o n red shift evolution found by lSerjeant et al.l (l2010l ) 



where SFR oc (1.0+z)^'^. ISerjeant et al.l (120101 ) did find that the highest luminosity quasars 



had a much stronger redshift evolution, perhaps as strong as SFR oc (1.0+2;)^*^. For simplicity 
the total volume-averaged star formation rate redshift evolution is used here, independent of 
the object luminosity. The local star formation rate is set such that SFR(z=0.0) = 0.5 
yr~^, in keeping with the average type 1 SFR found in a sample of local SDSS quasars by 



Kim et al.l (120061 ). Following the unified scheme, the SFR is not dependent on spectroscopic 



type. 

As shown by the red lines in figures El [8], [H [TOl [TTl and [121 the Euclidean normalized 
differential number counts increase with increasing flux level until peaking at ~1 mJy. The 
number counts decrease at flux levels above ~1 mJy, and, for wavelengths greater than 
160 /im, the number counts decrease more steeply as the wavelength increases. CT AGN 
dominate the counts on the brighter side of the peak and type 2 AGN dominate the lower 
flux level side of the peak. 

The star formation redshift evolution model contribution to the CIRB is shown in red 
in figure [13 The AGN contribution to the CIRB at 160pm is ~5%. CT AGN dominate the 
AGN contribution to the CIRB at wavelengths greater than ~100 pm, with type 2 AGN 
dominating at lower wavelengths. The AGN contribution to the CIRB for this model peaks 
at ~320 pm, a significantly longer wavelength than the peak of the CIRB as a whole. 

The P/^^*^ and p^iji^'^ for the star formation redshift evolution model are very similar 
to those for the constant star formation scenario, as shown by the red lines in figure [TH The 
total p/R of the star formation redshift evolution model under-predicts the observed local 
p/R but evolves more strongly with redshift than the constant star formation scenario. 
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4.4. Evolution with Redshift and AGN Lx 

The final star formation scenari o investigated here is the redshift and AGN Lx depen- 



dent star formation scenario used by IWilman et al.l (l2010l ). where 



SFR oc V^x/1043(i.o + 2)1-6. (7) 
Type 1 AGN are given the normalization constant .63 Mp^ yr~^ and fo r type 2 AGN the 



normalization prefactor is increased to 2.0 Mq jt~^ (jWilman et al.ll2010l ). High Eddington 
ratio CT AGN are given the type 2 SFR and low Eddington ratio CT AGN are given the 
type 1 SFR, for the reasons discussed in §4.21 

The Euclidean normalized differential number counts for the redshift and AGN Lx 
dependent SFR model are shown as blue lines in figures [71 El lOLfTOl fTTl and [I2] and are in 



EL HOI 
(120101 



decent agreement with the predictions made by IWilman et al.l (12010) based on a simulation 
of the extragalactic radio sky. The differential number counts increase with increasing flux 
level until peaking at 1-3 mJy, depending on wavelength. The peak of the number counts 
appears to increase with wavelength, with the peak at ~1 mJy at 100 fim and ~3 mJy at 
350 /xm. The number counts decrease from the peak to the brighter flux levels. At 100, 350, 
and 500 /zm CT AGN dominate the number counts on the brighter side of the peak, but at 
70, 160, and 250 /im type 2 AGN dominate except for at the brightest flux levels, >10 mJy. 

The blue lines in flgure [131 show the AGN contribution to the CIRB for the redshift and 
AGN Lx dependent SFR model. This model predicts ~5% of the CIRB at 160 /xm is due to 
AGN. Type 2 AGN dominate the AGN contribution to the CIRB at all wavelengths and CT 
AGN make a signiflcant contribution. For this model the AGN contribution to the CIRB 
peaks at ~300 fim, which is signiflcantly different from the peak of the total CIRB around 
160 /im. 

In flgure HH the pm for the redshift and Lx dependent SFR model is shown in blue. 
This model is in decent agreement with the pYr^'^ cit all redshifts but over-predicts the local 
P^R^'^ and under-predicts the local total piR. Both the total piR and p^R^^^ predictions 
are in agreement with observations at moderate redshifts, but do not evolve strongly enough 
with redshift to be in agreement with observations at the highest redshift. 
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Discussion 



We have presented predictions for observat ions of AGN and AGN hos ts in the FIR 
Herschel bands based on the composite model by Draper fc Ballantyne ( 2010 1^ These find- 
ings demonstrate that while AGN may not contribute a large fraction of the CIRB, AGN 
will be significant FIR sources and care must be taken in FIR surveys to identify AGNs 
as such. Here we discuss the implications of these results in terms of AGN and AGN host 
demographics. 



5.1. CT AGN 



A substantial population of CT AGN are necessa ry for AGN population synthesis models 



to ma t ch the peak of the CXRB at ~30 keV (e.g., iBallantyne et al.ll2006al : iTreister et al. 
2009al : iDraper fc Ballantynd 120101 ). Due to the extreme levels of obscuration in CT AGN, 
these elusive sources are generally only observed in the very hard X-ray, > 10 keV, or 
the infrared, especially the FIR. Using X-ray stacking methods, it has been shown that a 
large fraction of bright i nfrared excess sources {f^Alfn > 1000) host heavily obscured AGN 



e.g.. iDaddi et al.l 120071 : iFiore et al.l 120091 : ITreister et al.l l2009b[ ). However, it is uncertain 



how many of these highly o bscured AGN are actually CT (see iGeorgakakis et al.l 12010 



Georgantopolous et al.ll2010l ). Also, there is much de bate over the pr e valence of AGN in 



sources with more moderate infrared l umin osities (e.g.. iDey et al.ll2008l : iDonley et al.ll2008 



Fiore et al.ll2009l : ITreister et al.ll2009bll2010( ). Because CT AGN are generally not observable 
in the 2-10 keV band, AGN hard X-ray luminosity functions do not include CT AGN and 
therefore they must added in by hand to population synthesis models. In the population 
synthesis model used in this study, CT AGN are assumed to be accreting at either at L/Leh a 
> 0.9 or L/LEdd < 0.01, with fcr independent of /2, as found by lDraper fc Ballantynd ( l2010l ). 
Here we discuss predictions specifically for CT AGN in the FIR. 

At all wavelengths and for all star formation scenarios, the differential number counts 
are dominated by CT AGN for fluxes > 1-10 mJy. Type 2 AGN dominate at lower fluxes. 
Moreover, CT AGN, should constitute a non-trivial fraction of Herschel sources. Depending 
on the star formation scenario considered, CT AGN could make up ~10% of Herschel sources, 
even at 500 /im. 

Depending on the star formation scenario, CT AGN are found to contribute <5% of the 
CIRB at 160 /im. For bare AGN, CT AGN dominate the AGN contribution to the CIRB 



'^Models are available in table form by contacting the authors. 
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at wavelengths > 200 fim. When star formation is included, CT AGN dominate the AGN 
con tribution to th e CIR B at wavelengths > 100 fim for all star formation scenarios, except 
the IWilman et al.l ( 120101 ) model. 



Comparing the pr edictions made here against the infrared luminosity density found by 
Le Floc'h et al.l (120051 ). we find that CT AGN and their hosts contribute ~3% of the local 
infrared luminosity density from sources with Lj^ > 10^° Lq and nearly one-fourth of the 
infrared luminosity density from sources in the ULIRG range. In the AGN evolution star 
formation scenario, CT AGN can account for all of the local ULIRG range luminosity density. 
At 2; ~ 1 the relative CT AGN contribution decreases significantly in all luminosity ranges. 
However, when t aking into acc o unt t he stronger evolution of the SFR in high luminosity 
sources found by lSerjeant et al.l (120101 ). CT AGN and their hosts can still contribute nearly 
a quarter of the infrared luminosity density in the ULIRG range at 2; ~ 1. Showing that at 
higher redshifts CT AGN contribute less to the total infrared luminosity density, but may 
still contribute quite significantly to the brightest sources. 

The FIR is an important wavelength range for observing CT AGN due to the large 
amount of cold dust which obscures CT AGN. The majority of AGN observed by Her- 
schel will be CT. Depending on the star formation trends in CT AGN hosts, CT AGN and 
their hosts may constitute nearly ~10% of Herschel sources at 500 fim. The relative con- 
tribution of CT AGN and their hosts to the ULIRG range infrared luminosity density is 
5,25% and appears to be ap proximately constant over the redshift range z = 0-1 . However, 
Hatziminaoglou et al.l ( 12010! ) showed that AGN cannot be identified by their Herschel-SPIRE 
colors alone. Therefore fi nding CT AGN in the FIR Av ill require either Spitzer-MlPS coverage 
of bright SPIRE sources (IHatziminaoglou et al.ll2010l ) or X-ray stacking. Since the AGN and 
host differential number counts for both the composite and original models are dominated 
by CT AGN in the SPIRE bands. X-ray stacking of bright SPIRE sources is likely to disclose 
a large fraction of the CT AGN population. 



5.2. Differences Betvi^een Original and Composite Model 

The difference between the original model and the composite model is that in the com- 
posite model the CT AGN are put in specific, physically motivated Eddington ratio bins. 
Also, in the original model fcx /2, but in the composite model fcr is independent of f2- 
In order to understand the effects of the differences between the two models, we compare the 
predictions of the original and composite model for the constant star formation scenario. 

In the differential number counts, the differences between the composite and original 
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models are small but not insignificant. At all wavelengths the original model has a steeper 
decline in the bright end counts than the composite model. The original model predicts 
that the number counts will be dominated by type 2 AGN except for at the brightest fiuxes. 
Conversely, in the composite model the differential counts are dominated by CT AGN at 
every wavelength for fluxes > 1 mJy. 

The original model predicts a smaller overall AGN contribution to the CIRB than the 
composite model. In the composite model, CT AGN dominate the AGN contribution to 
the CIRB at wavelengths greater than 100 /im, but in the original model the type 2 AGN 
dominate the AGN contribution to the CIRB at all wavelengths. 

The original and composite models make similar predictions as to the AGN contribution 
to the infrared luminosity density in all luminosity ranges. The difference between these 
two models is most noticeable in the contribution of CT AGN to the local ULIRG range 
luminosity density. The original model predicts that CT AGN contribute 4% of the local 
ULIRG infrared luminosity density. The composite model predicts that nearly one-fourth of 
the local ULIRG infrared luminosity density is due to CT AGN. 

The overall predictions of the original and composite models are in agreement. To 
observationally determine which model best describes the CT AGN population, rigorous 
measurements of the CT AGN contribution to the local infrared luminosity density and/or 
accurate accounting of Herschel number counts to see whether CT AGN or type 2 AGN 
dominate the number counts will be needed. 



5.3. Star Formation in AGN hosts 

As seen in flgures [71 [8], [9l [TOl [11], and [12], the star formation scenario used greatly affects 
the predicted differential number counts. At wavelengths < 250 fim, the faint end slope is 
similar for all the star formation models investigated here, but the bright end slope is highly 
dependent on the host star formation at all wavelengths. The differences between the various 
star formation scenarios becomes more prominent when observing at longer wavelengths. The 
constant star formation model peaks at ~ 1 mJy and has a relatively flat bright end slope. 
The AGN evolution star formation model peaks at > 10 mJy with a knee at ~ 1 mJy. The 
redshift only evolution model peaks around, or just short of 1 mJy. This model also has the 
steepest faint and b right end slopes of the star formation scenarios considered here. The 
Wilman et al.l (l2010f ) star formation model peaks between 1 and 10 mJy for wavelengths > 



100 nm. The flux level of the peak in the differential counts may be an important tool in 
understanding the evolution of SFR in AGN hosts. This tool will be most effective when 
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observing at longer wavelengths. 



Dust obscured star formation is believed to be the primary progenitor of the CIRB 
with much debate as to the contribution from AGN. iBallantyne &: PapovichI (120071 ) found 
that AGN and star formation in AGN hosts can account for ~30% of the CIRB at 70 



/im. Similarly, Mullanev et all (l2010f l find that AGN contribute 5-25% of the CIRB at 70 



(Jauzac et al. 


, — ^ 
2010; 


Lacev et al. 


2010) 



formation contribute ~5% of the CIRB at 160 nm. Howeve r, when investigating the submil 



limeter properties of X-ray-selected AGN, iLutz et al.l ( 120101 ) found the average AGN SFR to 
be ~30 Mq yr~^, in which case AGN and their hosts would contribute ~88% of the CIRB 
at 160 yum. Therefore, understanding the star formation trends in AGN hosts is necessary 
for understanding how significant the contribution of AGN and their hosts is. 

The peak intensity of the CIRB occurs at ~160 fim. The star formation scenarios 
investigated here which take into account the evolution of the SFR with redshift predict the 
peak of the AGN contribution to the CIRB occurs at ~300/im. This suggests that the SFR 
of AGN hosts evolves differently with redshift than the SFR of normal galaxies. This effect 
could also be explained if the star formation SEDs used here are on average too cold, which is 
unlikely. If it is true that the SFR of AGN hosts evolves differently than the SFR of normal 
galaxies, this could offer important insights into the role of the AGN in the host galaxy 
evolution. At wavelengths shorter than the peak of the CIRB, the dominate contribution is 
from sources z < 1. Sources at z > 2 tend to dominate at wavelengths longer than 500 /xm. 

Another tool used to study differ ent AGN host star fo rmation scenarios is the infrared 
luminosity density. Using simulations, iHopkins et al.l ( 120101 ) find that AGN contribute 1-5% 
of the total i nfrare d luminosity density at all redsh i fts. Ap plying the classification scheme of 
Yuan et all koid l to AKARI sources, ICoto et all J2010bh find that AGN contribute ~20% 
of the total infrared luminosity density, ~40% of the luminosity density from sources in the 
LIRG range, and >90% of the luminosity density of sources in the ULIRG r ange, regardless 



of red shift. When considering the infrared luminosity density as measured by lLe Floc'h et al. 



( 120051 ). we find that the AGN and host galaxy contribution to the local infrare d lumin osity 
density is approximately a factor of 2 smaller than that found by iGoto et al.l ( l2010bl ). At 
2 ~ 1 the relative AGN contribution decreases to <5% for all luminosity ranges. Even when 
taking into ac c ount t he stronger evolution of the SFR in high luminosity sources found by 
Serjeant et al.l (120101 ) . it appears that AGN and their hosts only contribute ~25% of the 
infrared luminosity density in the ULIRG luminosity range at z ^ 1. The reduction in the 
AGN contribution to the ULIRG range infrared lumi nosity density by a factor of 2 between 
z ^ and 2; ~ 1, is consistent with the findings of ISturm et al.l (|2010[ ) that ULIRG level 
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luminosities can be achieved without major mergers at higher redshifts, suggesting that the 
AGN fraction in the high redshift ULIRG population will be smaller than that found locally. 

Determining the star formation history of AGN hosts is important in understanding 
why some galaxies host active supermassive black holes (SMBHs) and other galaxies host 
inactive SMBHs. The predictions presented here show that the flux level of the differential 
number counts peak in longer wavelength bands will be a helpful tool in determining the 
star formation history of AGN hosts. By comparing the peak wavelength of the AGN 
contribution to the CIRB to the peak wavelength of the CIRB, it is possible to determine if 
the star formation evolution of AGN hosts is different from that of normal galaxies. Also, 
at higher redshift AGN will have a smaller contribution to the ULIRG population, and that 
contribution will be dominated by CT AGN. 



5.4. Eddington Ratio Breakdown 

By using d^\(Lx , z) /d loRLx, the evolving Eddington ratio space density computed by 



Draper &: Ballantynd (|2010| ). instead of a traditional luminosity function, it is possible to 
make predictions for the contribution of AGN with different Eddington ratios. In figures 
HEl and dni the composite model Euclidean normalized differential counts are shown for the 
AGN evolution star formation model at 160 and 500 fim with the relative contributions from 
the different Eddington ratio bins shown. The blue lines show the contribution from AGN 
with L/LEdd < 0.01. The green lines show the contribution from AGN with 0.01 < L/LEdd 
< 0.9. The contribution from high Eddington ratio sources is shown in red. For all but the 
AGN evolution scenario, the number counts are dominated by low Eddington ratio AGN at 
all flux levels, due to the high space density of low accretion rate AGN at all redshifts. As 
shown in figures [T5] and [T6| the AGN evolutionary scenario bright end counts are dominated 
by high Eddington ratio CT AGN because of the high star formation rate in these objects. 

For all star formation scenarios the low Eddington ratio sources dominate the AGN 
contribution to the CIRB at wavelengths > 100 /im. The AGN contribution to the infrared 
luminosity density is dominated by moderate Eddington ratio AGN at all luminosity levels 
for all star formation scenarios except the AGN evolution star formation model, where the 
high Eddington ratio sources dominate the ULIRG range. The star formation scenarios that 
do not include redshift evolution find that low Eddington ratio sources dominate the AGN 
contribution to the local infrared luminosity density. 
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5.5. Implications for Herschel and ALMA 

The three wide area surveys conducted by Herschel will yield a large catalog of AGN host 
galaxies. Using the planned survey depths and areas as described in Section [H the models 
discussed here predict the following numbers of AGN hosts to be observed by the Herschel 
wide field surveys. Not including the planned lensing cluster observations, PEP should 
observe 100-500 AGN hosts, depending on the star formation scenario, at 160 fim with 5a 
significance. The portion of HerMES conducted during the Herschel science demonstration 
phase should provide 140-2100 AGN hosts, depending on the star formation scenario, with 
5cr significance at 250 fim. The H- ATLAS survey will yield 250-4000 AGN hosts at 160 fim 
and 90-1200 AGN hosts at 350 fxra with 5a significance. The H-ATLAS observed AGN host 
counts drop to 30-230 in the 500 /im band. This catalog of sources will provide a robust 
sample to better constrain the star formation properties of AGN hosts. 



Hatziminaoglou et al.l (120101 ) showed that AGN cannot be differentiated from normal 
galaxies based only on their FIR colors; therefore, ALMA will only be able to offer sup- 
plementary data on AGN which are identified in other wavelength bands. The wide area 
surveys conducted by Herschel will provide many promising candidate targets for ALMA. 
The dotted grey lines in figures [11] and [12] show the sensitivity limit of the 350 and 450 fim 
ALMA bands, respectively, as quoted in section [TJ Based on this projected sensitivity and 
the models presented here, the areal density of AGN hosts available to ALMA at 450 fim 
will be 300-1500 deg~^, depending on the star formation scenario, for an integration time 
of only 60 seconds. Therefore, using ALMA to conduct follow up observations on X-ray 
selected AGN should be an efficient way to study star formation in AGN host galaxies. As 
discussed in section 15. 3^ the differential number counts for different star formation scenarios 
seem to peak at different flux levels. These peak fluxes become more differentiated at higher 
wavelengths. By providing deeper submillimeter observations of Herschel sources and X-ray 
selected AGN, ALMA should be able to determine the evolution of star formation in AGN 
hosts. This will allow the determination of the wavelength where the AGN contribution to 
the CIRB peaks, and therefore whether the star formation in AGN hosts evolves differently 
than in normal galaxies. Furthermore, since ALMA will be taking spectra, measurements 
such as gas velocity, abundances, and temperature can be made. This will allow ALMA 
to not only determine the evolution of star formation in AGN hosts, but also to probe the 
physical structure of AGN host galaxy star formation. 
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5.6. Evolution of AGN 



From deep X-ray surveys it is known that high luminosi ty quasars and na oderate lumi- 



nosity AGN evolve differently with respect to redshift (e.g., lUeda et al.ll2003l ). This would 
suggest that quasars and Seyferts are caused by processes with different time scales. Also, 
it appears that the Seyfert population is well described by the unifi ed scheme, but the high 



luminosity AGN population may not follow the unified model (e.g.. iBallantyne et al.ll2006b 



Lutz et al.ll2010l ). A picture is starting to surface where the most powe r 



gered by major mergers as explored in the simulation s by 



ul A GN are trig- 



Hopkins et al.l (120061) and others 



but most AGN are triggered by less violent processes (IBallantvne et al.ll2006a[). Th i s pict ure 
is taken into account here by using the composite model of iDraper fc Ballantynd ( 120101 ) to 
describe the CT AGN fraction as Eddington ratio dependent. The difference in triggering 
process will also affect the star formation within the host galaxy. Mergers will engender not 
only AGN activity but also bursts of star formation. Where as the more secular evolution 
experienced by moderate power AGN will induce moderate, constant star formation. 

This AGN evolution scenario was further explored here with the AGN evolution star 
formation scenario. In this scenario, high Eddington ratio CT AGN hosts were assigned SFR 
= 175 M0 yr^^. The differential counts for this star formation scenario show a large peak 
> 10 mJy. However, the high Eddington CT A GN account for > 100% of the local ULIRG 
range infrared luminosity density measured by iLe Floc'h et al.l (j2005[ ). If the SFR in high 
Eddington ratio CT AGN hosts is reduced to 100 M© yr'^ then AGN will contribute 80% 
of the local ULIRG range infrared luminosity density, with a larger contribution to the local 
LIRG range infrared luminosity density. The prominent peak in the differential counts at ~ 
10 mJy persists, despite the reduced SFR in high Eddington ratio CT AGN. Therefore, a 
test of the AGN evolution scenario will be if the observed differential counts feature a strong 
bright end peak, possibly with a knee at moderate fluxes. As a large fraction of the sources 
contributing to this peak will be CT AGN, the AGN samples used for this test will have to 
be chosen very carefully. 



6. Summary 

FIR observations by telescopes like Herschel and ALMA will provide important insights 
on many questions about AGN and their hosts. Determining the flux level at which the 
differential AGN and host number counts peak will offer crucial constraints to the star for- 
mation history of AGN hosts, especially when observing at wavelengths at ~500 fim. The 
predictions presented here show that it is likely that the SFR in AGN hosts evolves differ- 
ently than the SFR in normal galaxies, as indicated by the peak wavelength of the AGN 
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contribution to the CIRB being significantly longer than the peak wavelength of the CIRB. 
Understanding how the SFR evolution in active galaxies differs from quiescent galaxies will 
provide clues on the triggering mechanisms of AGN and how the AGN interacts with the 
host galaxy. FIR observations will also allow the AGN evolution scenario to be tested by 
comparing SFRs in bright AGN with different levels of obscuration. The relative contribu- 
tions of AGN with various levels of obscuration to the bright end differential counts will also 
be an important test of the major merger trigger model. Applying X-ray stacking techniques 
to bright 350 or 500 //m sources, especially sources with a hot dust component in the SED, 
will be an efficient way of finding CT AGN. 

The authors thank A. Poglitsch and I. Valtchanov for information on the PACS and 
SPIRE filters, P. van Hoof for assistance with Cloudy, and the referee for a useful report 
that helped to improve this paper. This work was supported in part by NSF award AST 
1008067. 
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Fig. 1.— Rest frame SEDs for high Eddington ratio AGN {L/LEdd > 0.9) with logL^ = 43 
and z = 0.45 (fcT = 0.86 and /2 = 0.78). The type 1 SED is shown in red, the type 2 SED 
is shown in green, while the Compton thick SED is shown in blue. 
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Fig. 2. — Lx versus L^.s^im for the high Eddington ratio {L/LEdd > 0.9) model SEDs at 
z = 0.05. Type 1 SEDs are shown as stars, type 2 SEDs are shown as diamonds, and 
the Compton thick SE Ds are shown as tria ngles. The black line is the best-fit line for the 
well- resolved sample of iGandhi et al.l ( l2009l ). 
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Fig. 3. — Euclidean normalized differential number counts for bare AGN at 70, 100, 
160, and 250 /im. The bla ck lines plot the predictions based on the composite model of 



Draper &: Ballantynd (l2010l ) and the cyan lines show the predictions for the original model. 
The dotted-lines are predictions for type 1 AGN, the dot-dashed lines are predictions for 
type 2 AGNs, and the dashed lines are predictions for Compton thick AGN. 
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Fig. 4. — Comparison of Euclidean normalized differential number counts for bare AGN at 
70, 100, 160, and 250 /im using the composite model. The blue lines are r = 1 pc and uh ct 
= 10^ cm~^, the black lines are r = 10 pc and uhct = 10^ cm~^, and the green lines are 
r= 10 pc and nH,cT = 10^ cm~'^. The line styles are the same as in figure [3l 
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Fig. 5. — Bare AGN contribution to the CIRB. Line color and styles are the same as 
in figure [HI Data points are from a variety of instruments: the fill ed diamonds are from 
lACTS jMazin fc BauA boOTl ): the open diamoiid is f rom MIPS JPapovich etaP UQM)', 
the triangles are from DIRBE Jpinkbeiner et all I2OO0I): the pentag ons are from AKARI 
(IMatsuura et al.ll2010l): the asterisks are from MIPS (IDole et al.ll2006r): the s quares are from 



DIRBE flHauser et al. 



1998r): the cross is from ISOPHOT (jjuvela et al.ll2009l): the star is also 



from ISOPHOT (iLagache fc Pugetll2000l ): and the circles are from BLAST (jPatanchon et al. 
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Fig. 6. — Infrared luminosity density of bare AGN with respect to redshift. Sohd hnes and 
circles refer to pm, dot-dashed lines and triangles refer to pYr^*^ ■, while dashed lines and 
squares refer to p^li^^^ ■ The black lines show the composite m odel and the cyan li nes show 
the original model. The z = 0.0082 data points ar e taken from lGoto et al.l (j2010bl ) and the 
higher redshift data points are from the work of iGoto et al.l (j2010al ) in the AKARI NEP 
deep field. 
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Fig. 7. — Euclidean normalized differential number counts for AGN and host star formation 
at 70 /im for the various star formation scenarios using the composite model. Constant 
star formation is shown as black. The AGN evolution st ar formation sc e nario is shown as 
green. Star formation with the redshift evolution found by ISerjeant et al.l (120 lOf) is shown as 



red. T he star formation scenario using the redshift and AGN Lx evolution by IWilman et al. 



( I2OIOI ) is shown in blue. The line styles are the sam e as in figure [3l Th e grey line shows the 
best fit model galaxy differential numb er counts of iFranceschini et al.l ( 120101 ). Spitzer data 
points are from iBethermin et al.l ( I2OIOI ). 
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Fig. 8. — Euclidean normalized differential number counts for AGN and host star formation 
at 100 fim for the various star formation scenario s using the compos ite model. Lines are the 



same as in figure [71 Herschel data points are from lAltieri et al.l (120101 ) (stars) and lBerta et al. 
fcoioh (circles-GOODS-N, squares-Lockman XMM, and triangles-COSMOS). 
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Fig. 9. — Euclidean normalized differential number counts for AGN and host star formation 
at 160 fim for the various star formation scenarios using the composite model. Line colors 
and styles are the same as in figure [71 Data points are the same as in figure El 
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Fig. 10. — Euclidean normalized differential number counts for AGN and host star formation 
at 250 fim for the various star formation scenarios using the composite model. Line colors 



and s tyles are the same as in figure [71 Circles are data points f rom Hersche l (lOliyer et al. 



2010l ) and astrisks show the multiply-broken power-law model of iGlenn et al.l (120101 ). 
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Fig. 11. — Euclidean normalized differential number counts for AGN and host star formation 
at 350 fim for the various star formation scenarios using the composite model. Line colors 
and styles are the same as in figure [71 Additionally, the dotted grey line shows the expected 
continuum sensitivity of ALMA for an integration ti me of 60 seconds an d a spectral resolution 
of 1 km/s. Circles are data p oints from Herschel ( jOliver et al.ll2010l ). the diamond is from 
SHARC II (IKhan et al.l 120071 ). and astrisks show the multiply-broken power- law model of 



Glenn et all tOld) . 
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Fig. 12. — Euclidean normalized differential number counts for AGN and host star formation 
at 500 fim for the various star formation scenarios using the composite model. Line colors 



and s tyles are the same as in figure [TTl Circles are data points fro m Herschel { Oliver et al. 



2010l ) and astrisks show the multiply-broken power-law model of I Glenn et al.l (120101 ). At 



bright fluxes numerical artifacts are present due to the small number of sources in this flux 
region. 
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Fig. 13. — AGN and host star formation contribution to the CIRB for the various star 
formation scenarios using the composite modeL Line colors are the same as in figure [HI The 
hne styles are the same as in figure El Data points are the same as figure [5l 
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Fig. 14. — AGN and host star formation infrared luminosity density for the various star 
formation scenarios using the composite model. Line colors are the same as in figure [HI The 
line styles and data points are the same as figure El 
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Fig. 15. — Euclidean normalized differential number counts for AGN and host star formation 
at 160 fim for the AGN evolution star formation model. Black lines show the total predictions 
for the composite model with the low Eddington ratio sources in blue, the mid Eddington 
ratio sources in green, and the high Eddington ratio sources in red. The line styles are the 
same as i n figure [3l The solid grey lines shows the best fit model galaxy differential number 
counts of iFranceschini et al.l (120101 ). the Data points are the same as in figured 
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Fig. 16. — Euclidean normalized differential number counts for AGN and host star formation 
at 160 fim for the AGN evolution star formation model. Line colors and styles are the same 
as in figure [15] with the addition of the dotted grey line which shows the expected continuum 
sensitivity of ALMA for an integration time of 60 seconds and a spectral resolution of 1 
km/s. Data points are the same as in figure [TOl At bright fluxes numerical artifacts are 
present due to the small number of sources in this flux region. 



